Quantum information approach to the quantum phase transition in the Kitaev 

honeycomb model 
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Kitaev honeycomb model with topological phase transition at zero temperature is studied using 
quantum information method. Based on the exact solution of the ground state, the mutual infor- 
mation between two nearest sites and between two bonds with longest distance are obtained. It is 
found that the mutual information shows some singularities at the critical point the system transits 
from gapless phase to gapped phase. The finite-size effects and scaling behavior are also studied. 
Our results indicate that the mutual information can serve as good indicator of the topological 
phase transition. This is because that the mutual information is believed to be able to catch some 
global correlation properties of the system. Meanwhile, this method has advantages that the phase 
transition can be determined easily and the order parameters, which are hard to be obtained for 
some topological phase transitions, are not necessarily known. 

PACS numbers: 03.67.-a 05.70.Jk 05.30.Pr 75.10.Jm 



I. Introduction 

Recently, Kitaev honeycomb model has become a pop- 
ular subject in both the fields of condensed matter 
physics and quantum information processing [ll-fl3j . This 
model was first introduced by Kitaev to study the anyons, 
and the analytic exact solution to the ground state of 
this model has been obtained by several methods [l|, \\a\ - 
\vf\ . It has rich phase transitions and has both a gap- 
less phase with non-Abelian anyons excitation and three 
gapped phases with Abelian anyons excitations depend- 
ing on the values of the parameters in the Hamiltonian. 
It is shown that the system possess a topological phase 
transition, which is not able to be characterized by sym- 
metry breaking theory and the corresponding local order 
parameters but can be characterized by nonlocal string 
order parameters p], Il8| . One interesting point is that 
the system is a scarce exactly solvable model with di- 
mensions higher than one , thus it provides a test bed for 
many numerical methods in two dimensional systems just 
as the Ising model does in one dimension. With these in- 
teresting properties, Kitaev honeycomb model has been 
studied intensively and is extended to other cases [l^.ll9j. 

On the other hand, Kitaev honeycomb model also has 
good practical advantages to be an active subject in that 
it has great potential applications in the quantum in- 
formation and quantum computation. It was suggested 
to use the Kitaev honeycomb model to realize the fault- 
tolerant topological quantum computation. The system 
is a good candidate to encode quantum information while 
those quantum states can be naturally protected from 
the inevitable decoherence by environment poj . The Ki- 
taev honeycomb model can be realized using the optical 
lattice [2ll |22|. and using the superconducting quantum 
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circuits [23|, |24|. It is also studied by means of fidelity 
susceptibility (1 5| and the extended Kitaev model is stud- 
ied by approaches of entanglement [l3| . 

In this paper, we investigate the Kitaev honeycomb 
model from the quantum information perspective |26| . 
We study the topological phase transition in this model 
by means of mutual information between the component 
lattices. It is generally believed that the mutual infor- 
mation measures the total information and describes the 
global correlation properties [2 71]. We find that both the 
derivative of mutual information between two nearest 
neighbor lattices and the mutual information between 
two bonds of the lattice can detect the topological phase 
transition in the Kitaev honeycomb model. This quan- 
tum information method has great advantages in that 
the singular behavior occurs exactly at the point when 
the gapless phase transits into a gapped phase. We also 
study the finite-size effects and the scaling behavior of 
the singularities of the mutual information. 

This paper is organized as follows. In Section II we 
briefly introduce the Kitaev honeycomb model, then di- 
agonalize the Hamiltonian and give the exact solution of 
the ground state based on the initial Kitaev's method. 
After that, we calculate the two sites and four sites cor- 
relation functions getting prepared for the two-site and 
two-bond reduced density matrix. In section III and sec- 
tion IV, we calculate the two-site mutual information 
and two-bond mutual information and the former one's 
derivative, respectively. Section V is the conclusions and 
remarks. 



II. Kitaev honeycomb model 



Kitaev honeycomb model is a two dimensional spin- ^ 
lattice model with nearest neighbor interactions. It has 
two kinds of simple sublattice which are denoted by the 
dark dots and empty circles in figure 1. Each lattice in- 
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FIG. 1: (color online). The sketch map of Kitaev honeycomb 
model. The unit cell contains two sites of different kinds, 
which is highlighted by a elliptic circle. For simplicity, we 
choose the coordinate axes in ni and ri2 directions. 



teracts with three nearest neighbors of the opposite kind 
through three distinct bonds labeled as x link, y link 
and z link. For each bond the interaction has different 
coupling constant. The Hamiltonian is 



H 



-Jx E tf^k-Jy 



x — links 



E 

y— links 



v v 



J - E a 



-links 



(i) 

where the subindex j, k denote the location of the site, 
and cr^(a = x,y,z) is the pauli matrix at site k. We 
take axis of the system in the and n 2 direction, and 
in each direction there are L unit cells. Therefore the 
whole system has 2L 2 sites. Next, we used the original 
Kitaev's method to diagonalize this Hamiltonian and get 
its ground state. 



A ^ B. Thus, the Hamiltonian becomes 

H = -E JQ E b PkW 

ot ex.— links 

= *E JQ E Wte<* 



a a— links 

*E JQ E ujkCjCk 

ot a— links 



o E Ja i. 



k u hk c 3 c k . 



(3) 



fn the last equation, the value of a is totally determined 
by the site index j and k. The factor \ is due to the 
summation of the lattices has counted each lattice twice, 
ft can be easily shown that u 2 . k = 1, [uj : k,H] = 0, and 
Uj t k commute with each other. As a result the eigen- 
values of Uj,k here we present by Uj t k are ±1, and the 
whole Hilbert space can be decomposed into a series of 
eigenvalue spaces described by the eigenvalues of Uj t k- 
According to [13, HH , the ground state is in the vortex 
free space so that we assume Uj t k = 1 f° r ah links, where 
j is a kind of simple sublattice presented by the empty 
circles in this paper. Notice Uj k — —Uk,j- 

As the unit cell of this model contains one empty circle 
lattice and one dark dot lattice, we introduce a pair of 
index (s, A) to take the place of the previous site index j, 
where the first index s stands for the location of the unit 
cell, and the second one describes the two different kinds 
of sublattice. In this paper, we let the empty circle's 
second index takes the value 1, and the dark dot's takes 
the value 2. Please see figure 1. Then the Hamiltonian 
becomes 



S, A,t,/2 



(4) 



The two dimensional system we studied is on the surface 
of a torus with periodical boundary conditions. Because 
of the translational invariance of the system J s ,A,t,^ is 
actually determined by three index A, fi, and t — s. Then 
we introduce the Fourier transformation 

J r .,A,t, M = Jo,A ; t-s, M =^X) e ^ q ' (rt " ri) ^( t l)' 



A. The ground state 



C S ,A 



E 



e iqrs a q , A . 



(5) 



We first introduce the following Majorana transforma- 
tion to transform the Pauli operators into the Majorana 
fermion operators. 



The inverse transformation is 

J^(q) = 5> <q ' rtJ M,t,M. 



g x = ib x c, a v — ib v c, <r z — ib z c, 



(2) 



«q,A 



2L 2 E 



(6) 



where the Majorana operators satisfy — A, A 2 = 1, 
AB + BA = 0, b x b v b z c = I, for A,B E {b x , b y , b z ,c} and 



where 



a q ,A 

t 



satisfies 
,t 



-q,A 



^q,A' "q,A 



0, 



+ = a q,Aa qiAI + a q , M « q ,A 



S P q5\ ^, and other 



3 



anticommutators are all equal to zero. Then the Hamil- 
tonian arrives at 

2 

H = Y J A^(q) a -q,Aa q ,M- ( 7 ) 

After simple calculations we obtain that Ji,i(q) = 
J2 t e jq r * Joi.ti = 0, because Joi.ti = 0. For the simi- 
lar reason J 2) 2(q) = 0. As each lattice interacts with 
its three nearest neighbors, there are only three values 
of t corresponding to the three neighbors that make 
Jbi,t2 take nonzero values. Thus = J x e lcl ' ni + 

Jj,e iqn2 + J z , and J2,i(q) = —Ji,2 (q), where ni and n 2 
are in certain directions which is shown in figure 1. Let 
/(q) = Ji i2 (q) = e(q) + iA(q), and choose q x to be in 
the direction of ni, and to be in the direction of n 2 . 
Then we have 

e(q) = JxCosq x + JyCOsq y + J z , 
A(q) = J x sinq x + J y sinq y , (8) 

where q x and q y take values q Xl q y = 27rn/L, n = 
-(L-l)/2,--- ,(L-l)/2. We can see that e(-q) = e(q), 
A(-q) = -A(q), and /(-q) = /*(q). The Hamiltonian 
becomes 

# = Y i /( < 30 a q,i°q.2 + (*/(q))* a L fl q,i- ( 9 ) 
q 

Next, we introduce the following Bogoliubov transforma- 
tion: 

Cq.l = w q a q> i + v q a q , 2 , 

Ct * t r * t 

q,l = U q a q,l + U q<,2: 

Cq,2 = W q aq,l - u q aq,2, 

C l,2 = W q«q,l ~ W q°q,2' ( 10 ) 

with the new operators satisfying [Cq y \, „]+ = 
<5pq<^A,/i ,Cq _^ = 0. Using the Bogoliubov transformation, 
the Hamiltonian is diagonalized as 

H = Y1 \f<*K C li C «A -C q , 2 C q ,2), (11) 
q 

with Wq = Uq = -*=^, w _ q = - Uq , C_ qil = 
— 2u q u q C q2 . Considering the fact C ql C q) i = 1 — 
C_ q2 C_ q)2 , the Hamiltonian reads 

q 

= ^|/ q |(l-2(7 q)2 C q , 2 ). (12) 
q 

The normalized ground state is 

i G > = U C U°)> ( 13 ) 

q 

with C qj2 |0) = 0. The energy gap is 2min q {|/ q |}. 
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FIG. 2: (color online). The phase diagram of the Kitaev 
honeycomb model in the J x + J y + J z = 1 plane in the pa- 
rameter space. In the three shadow areas labeled by A x , A y 
and A z , the system is gapped with Abelian anyon excitation, 
and in the blank area labeled by B the system is gapless with 
non-Abelian excitation. In this paper, we focus on the red 
dash line J x = J y = (1 — J 2 )/2, where the critical point of 
topological phase transition is J z — 0.5. 



B. Phase diagram 

This ground state has two distinct phases in the pa- 
rameter space. In the region of \J X \ < \J y \ + \ J Z \, 
\J y \ < \J X \ + \Jz\ and \J Z \ < \J y \ + \J X \ it is gapless 
with non-Abelian excitation and in other regions it is 
gapped with Abelian anyon excitations [l4j. We focus on 
the J x + J y + J z = 1 plane. The phase diagram is shown 
in figure 2. In this paper, we investigate the behaviors 
of two-site mutual information and two-bond mutual in- 
formation in the phase transition from the gapless phase 
to a gapped phase along the red dash line in the phase 
diagram of figure 2. 



C. Correlation functions 

In this section we calculate the two-site and four-site 
correlation functions at the ground state of the systems 
which will be used to construct the reduced density ma- 
trix. Suppose the two nearest lattices to be studied arc 
linked by z-bonds. The correlation function between two 
nearest lattices is 

«i< 2 > = (^i^4Ee i(q+<,,) ' r ^iV,2) 
qq' 

= -^E e<(qW) ' r <<w,2>- 

q,q' 



4 



By using the relation 

( a qa<V,2> = (( U q C 'q.l+ t; q C q,2)(«q'C , q' 1 l-«q' C 'q',2)) 

Lx A. 

~ 2 q '- q '|/ q |' 



we obtain the correlation function 



«i<2> 



L \ " JfL = _L_ \ " 



L 2 



2i 2 



1 £ q 



/q + /-q 

l/ql 



(14) 



where £ q = |/ q | = ^ ' e\ + A q . 



The correlation function between two bonds as high- 
lighted by elliptic circles in figure 1 is 



where 



( CT r 1 ,l CT r 1 ,2 <7 r 2 ,l <7 r 2 ,2) = (^.l^^^^l^^Cr^lCn^CrjjC^^) 

4 £ e J ^ + ^ ) - ri e j(q 3+^)-r 2(ariiiarii2ar2iar ^ 2)j 



L 4 



qi,q 2 ,q 3 <q4 



(a ri 4 a ri,2ar 2 ,iar 2 ,2) - |y^( C '-q 1 ,2( C, -q 2 ,2 + Cq 2;2 ) (Cl^ - Cq 3 , 2 )Cq 4i2 ) 

1 /qi /q 3 



4|/ | |/ |(^2.-q3 5 qi,-q4 ^qi,-q 3 ^q2>-q4 ^qi,-q 2 ^q 3 >-q4)- 



Then, we arrive at 



((T z z z z \ Lf V A_ l (q 1 -q 3 )-(r 1 -r 2 ) _ \^ i( qi +q 2 )-(ri-r 2 ) _ \^ Ai /q 3 

= E /qi iXj : i/q 1 ,!" 3 (cos[(qi " qs) ' (ri " rz)] " 1} 



qi<q 3 



I A qi A q3 ^q 3 (cos[(qi _ q3) . (ri _ r2)] _ l) . 



L 4 



E 

qi.q3 



E„, E, 



qi ^q.3 



r 



(15) 



III. Mutual information between two neighbor 
lattices 



The reduced density matrix of two site i and j is 

Pi,i = lE 3 a,0= O (°>j)°>j, Where ^ = ^ = ^ 

a- 3 = a z and er° is identity. In the system (1), each site in- 
teracts with its three neighbors by different operators^, 
a y and cr z ), while each two linked sites together have only 
one kind of operator available, i.e., a a a a , with a corre- 
sponding to the type of their link. We find that only the 
correlation function along the link interacting direction is 
nonzero when we study the reduced matrix of two linked 
sites. That is to say if we consider the two lattices with 
z link (please see figure 1), all the correlations are zero 
except (a z a z ). Therefore, the reduced density matrix of 
this model has only diagonal elements, although the in- 



teractions of one site have three components. It indicates 
that the model (1) is much more like a classical system 
and similar to the Ising model. That may explain why 
this two dimensional model can be solved analytically. 



In order to show our results more clearly, we use the nu- 
merical method to diagonalize the Hamiltonian exactly. 
With the periodical boundary conditions, we diagonalize 
an eight lattices system which is the smallest subsystem 
on the surface of a torus, please see figure 3, and calculate 
all the 16 two-site correlation functions. Here we omit the 
subindex corresponding to the sites' position because the 
value is invariant due to the translational invariance of 
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this system. The explicit form of the Hamiltonian reads 
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FIG. 3: (color online). The smallest translation invariant sub- 
system of Kitaev honeycomb model with periodical boundary 
conditions. The eight interaction sites in the subsystem are 
highlighted by red numbers in the graph. The sites labeled by 
small black numbers are the repetitions of the 8 sites because 
of the periodical boundary condition and torus topology. 



—X i —X —X \ 



Hg — J x (<7f<73 + <76<7 4 T u 7 u t 

+ J z (a%a% + a x 4 a* + ala x 1 +a*a%). (16) 
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By using the periodical boundary conditions, the system- 
size of the system (|16|) can be extend to infinity. The 
main properties of the system are kept since all possible 
interactions are considered. The result is shown in figure 
4. We see that the correlation functions along the in- 
direction is non-zero. 
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FIG. 4: (color online). The ground state correlation functions 
of two-site with z-link in the Kitaev honeycomb model. We 
see that (<Tr,i°r,2) — h (^,1^,2) 7^ and others are zero. J z 
is in unit of J x + J v + J z . 



The explicit form of reduced density matrix of two sites 
with nearest neighbor is 



Pr,l;r,2 



/ i + «,i< )2 > 

i-«,i<,2> 

1 - (<J z rul a z ru2 ) 



V 











(17) 



l + «i<2>/ 



The eigenvalues of the reduced density matrix are Ai = 



A 2 = (1 - (cr z (T z ))/A and A3 = A 4 = (1 + (cr z (J z ))/4. Ev- 
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FIG. 5: (color online). Mutual information between two con- 
nected sites. The mutual information increases monotonically 
with the increasing of J z . However, there exist a certain value 
J™. If J < J™, the mutual information is a concave func- 
tion while if J > J™, the mutual information is a convex 
function. Thus, the first order derivative of the mutual infor- 
mation with respect to J z has a peak at the J™. J z is in unit 
of J x + J y + J z . 

ery eigenvalue corresponds the two-fold degenerate eigen- 
states. From the two sites reduced density matrix we 
can derive the reduced density matrix for one site as 7/2, 
where I is identity matrix. The mutual information of 
the two sites i and j is 

S(j;3) = 2 - 2ff (l^>)_ 2ff (l±i^) (18) 

where H(x) — — xlog 2 (a;). 

The the two-site mutual information along the line 
J x — Jy = (1 — Jz)/2 is shown in figure 5. We see that 
the mutual information increases monotonically with the 
increasing of J 2 . However, there exist a certain value 
J™. If J < J™, the mutual information is a concave 
function while if J > J" 1 , the mutual information is a 
convex function. Thus, the first order derivative of the 
mutual information with respect to J z has a peak at the 
J™. The derivative of the mutual information is shown 
in figure 6. From figure 6, we see that the first order 
derivative of the mutual information arrives at the max- 
imum value at the point J™ = 0.5, which exactly cor- 
responds the critical point J z = 0.5. We also find that 
the value of J™ is fixed when the system-size changes. 
The maximum value is a constant when the system size 
tends to infinity, as shown in figure 7. This is different 
from the Ising model where the second order derivative 
of entanglement entropy diverges at the critical point in 
the the thermodynamic limit [3l|, although their density 
matrixes are same in structure. 

From the quantum information perspective, the entan- 
glement measured by concurrence between two sites is 



4.5 




FIG. 6: (color online). The first order derivative of mutual 
information between two connected sites for the system-size 
L = 100. The derivative of the mutual information has 
a maximum exactly at the critical point J z = 0.5. The 
small peaks in the gapless region is due to that while cal- 
culating the (crr.i "?^) certain e q and Eq are both infinites- 
imal and some systemical error are included. The peaks 
get smaller when L increases. The subgraph shows that 
when the system-size tends to infinity, the maximum of the 
derivative of the mutual information tends to a constant as 
log 2 U^§^-) max ~ 3.87934 J = -0.00384L - 2.12174. J z is 
in unit of J x + J y + J z . 



zero since the density has only diagonal elements, while 
the entanglement between one site and all the rest sites is 

maximum entanglement, i.e., c = \J — Trpf) = 1, 

where c denotes the concurrence, pi is the reduced den- 
sity matrix of the particle in site i and d is the dimension 

of Pi 0,1131. 



IV. Mutual information between two bonds with 
longest distance 

In this section, we investigate the mutual information 
between two z linked bonds. Firstly, we need to calculate 
the density matrix. For arbitrary two z linked bonds at 
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FIG. 7: (color online). The mutual information between two 
z linked bonds with longest distance located on the torus sur- 
face for a given system-size. The small peaks in the gap- 
less phase is caused by the 0/0 error as before. The sig- 
nificant peak arises exactly at the critical point. The sub- 
graph shows that when the system-size tends to infinity, the 
peak value of the mutual information trends to a constant as 



log 2 (S(ri : r 2 ) 



-8.34363 xlO -13 ) 



-0.10637L- 18.8454. 



Here J z is in unit of J x + J v + J z 



Y\ and r 2 the density matrix is 

1 \ a a 
Pr u r 2 = 2^ a ri, 1^,2^2, 1 CT 





r2,2 



a, 0=0,3 

\y I ^r a ^r a ^ ^ \ 

X \ <J r 1 , 1^,2^2, l a r 2 ,2)- 



(19) 



The eigenvalues of this density matrix are (1 — 
A(a z a z a z a z ))/l&, (1 - 2(a z a z ) + (a z a z a z a z ))/16, (1 + 
2{a z a z ) + {a z a z a z a z ))/16. These eigenvalues correspond 
the 8 fold, 4 fold and 4 fold degenerate eigenstates, re- 
spectively. Then we obtain the mutual information be- 
tween two z linked bonds in ri and r 2 as 



Sin : r 2 ) = AH 

- 8H 

- ah 

- AH 



1 + {cr z a z ) 



AH 



1 - (a z a z ) 



1 - (a z a z a z a z ) 
16 

1 - 2{a z a z ) + (a z a z a z a z ) 
16 

l + 2(a z a z ) + (a z a z a z a z ) 
16 



(20) 



In order to reveal the long range correlation in the 
system, we study the bonds in the direction that contains 
the largest distance in the torus, which is marked by a line 
of elliptic circles in figure 1. Without losing generality, 



we choose the positions of the bonds as ri = (0, 0) and 
r 2 = (0.5L,0.5L) in the (ni,n 2 ) coordinate system so 
that the two bonds are longest separated. The mutual 
information between two z linked bonds in this direction 
is shown in figure 7. We find that the mutual information 
has a maximum at the critical point. This property is 
valid for different system-sizes. The peak values tend to 
a constant when the system-size L trends to infinity. 



V. Conclusions 

In this paper, based on the exact ground state of Ki- 
taev honeycomb model we have obtained both the re- 
duced density matrices of two nearest neighbor sites and 
between two z linked bonds with longest distance. We 
show that these density matrices have only diagonal el- 
ements, so that there is no entanglement between two 
local sites or two local bonds, but the nonlocal entan- 
glement between one site and the rest of the whole sys- 
tem is maximum. From quantum information theory, the 
ground state which is a multipartite state seems more like 
Greenberger-Horne-Zeilinger (GHZ) states other than 
W-type states. We have calculated the mutual informa- 
tion between two nearest neighbor sites and the mutual 
information between two z linked bonds with the longest 
distance in the torus topology. The first order deriva- 
tive of the former mutual information and the latter 
mutual information itself have peak at the point where 
the ground state transits from the gapless phase into a 
gapped phase. The above singular behavior serves as 
an exact and easily obtaining detector of the topolog- 
ical phase transition in the Kitaev honeycomb model. 
The so called localizable entanglement is related with the 
string order parameters and the hidden topological long 
range order in one dimension spin chains [3 2l |33|. More 
over, the topological phase and topological phase transi- 
tion have their roots in the hidden topological long range 
order and the string order parameters. Therefore, the 
investigation of the relation among the localizable entan- 
glement, topological phase transition and other related 
quantities in the exactly solvable two dimensional Ki- 
taev model may greatly enhance our understanding of the 
topological phase. The research in this direction should 
be further explored to extensively study the topological 
phase and topological order by means of quantum infor- 
mation. 
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